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ABSTRACT 


Theoretical  equations  for  unsteady,  compressible,  laminar  flow  in 
cylindrical  tubing  have  been  derived  for  pressure  measuring  systems 
comprised  of  up  to  three  tube  sizes  connected  to  constant  volume  trans¬ 
ducers.  The  basic  equations  have  been  solved  for  both  step  and  ramp 
pressure  input  forcing  functions.  Digital  computer  programs  have  been 
written  to  be  used  as  engineering  tools  for  rapid  determination  of  system 
characteristics.  Examples  of  the  ability  of  the  programs  to  compare 
tubing  systems,  optimize  diameters,  generate  time  responses  to  syn¬ 
thesized  forcing  functions,  and  determine  lag  times,  distortion,  and 
attenuation  are  presented  as  guides  for  prospective  users  of  the  programs. 
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NOMENCLATURE 


A  Internal  cross-sectional  area  of  tube,  ft^ 

a  Coefficient  in  Weber's  equation,  elsewhere  the  dimensionless 

condition  for  neglecting  acceleration  forces 

B  Dimensionless  geometric  ratio 

C  Dimensionless  geometric  ratio,  or  a  constant  when 

accompanied  by  subscript 

D  Dimensionless  geometric  ratio 
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d  Internal  diameter  of  tube,  ft 

e  Natural  base,  2.  7183 

F(p)  Function  of  pressure  only,  psf--*- 

f  Frequency,  cps 

f(t)  Pressure  forcing  function,  psf/sec 

g(t)  Transformation  function 

i  VT 


K 

Km,KT 

K* 

k 

SL 

m 

n 

P 

Po 

P 

R 

Re 

T 

Tl 

Tr 

t 

u 

V 

W(a,  x) 

W'(a,  x) 

w 

x 


Constants  as  defined 

Geometric  characteristics  of  tubing  system,  psf-sec 
Km  .  10-3 

Forcing  function  coefficient,  psf/ sec3 
Length  of  tube,  ft 

Forcing  function  coefficient,  psf/sec^ 

Forcing  function  coefficient,  psf/'sec 

Transducer  pressure,  psf 

Initial  transducer  pressure 

General  pressure,  psf 

Gas  constant,  1715  ft-lb/slug-°R 

Reynolds  number 

Temperature,  °F 

Temperature,  °R 

Reference  temperature,  524.4  °R 

Time,  sec 

Time  lag,  sec 

Known  solution  to  Riccati  equation,  psf--*- 

Transducer  volume,  ft3 

Weber  function 

d  W(a,  x) 
dx 

Mass  flow  rate,  slug/sec 

Distance  along  tube,  ft,  or  transformation  variable  as 
defined 
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y  Transformed  pressure  function 

T  Gamma  function 

S,Ap  Amplitude  of  pressure  change,  psf 

e  Error  in  final  pressure,  psf 

/u  Viscosity,  slug/ft-sec 

pr  Reference  viscosity,  3.  8  x  10-^  slug/ft-sec 

ir  3.  1416 

t  Period  of  forcing  function  sec, 


SUBSCRIPTS 

1 

2 

3 

12 

23 


avg 

comp 

f 

i 

min 
opt 
t,  orf 
trans 


First  tube  from  orifice 

Second  tube  from  orifice 

Third  tube  from  orifice 

Junction  of  first  and  second  tubes 

Junction  of  second  and  third  tubes 

Average 

Compromise 

Final 

Initial  or  as  otherwise  specified 

Minimum 

Optimum 

Orifice 

Transducer 
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Measurement  of  a  large  number  of  test  article  pressures  during  wind 
tunnel  testing  invariably  requires  that  the  pressure  transducers  be  located 
outside  of  the  test  cell  because  of  the  physical  dimensions  and  environ¬ 
mental  limitations  of  present  day  high  accuracy  pressure  transducers. 

The  result  is  that  the  long  lengths  of  tubing  that  are  required  to  connect 
model  orifices  to  the  transducers  act  as  delay  lines.  If  the  orifice  pres¬ 
sures  are  steady  values,  the  pressures  seen  by  the  transducers  will 
eventually  become  the  true  orifice  values,  but  only  after  an  interval  of 
prime  testing  time,  following  a  change  in  pressure  conditions,  has  been 
expended.  With  the  model  orifice  sizes  and  tubing  configurations  presently 
being  used  in  the  1 6  -  ft  wind  tunnels  (Propulsion  Wind  Tunnels,  Transonic 
(16T)  and  Supersonic  (16S)),  it  is  not  uncommon  for  from  10  to  20  percent 
of  the  prime  on-line  testing  time  available  for  a  test  program  to  be  devoted 
to  pressure  stabilization  time.  Following  a  Mach  number  extension  of  the 
supersonic  tunnel  to  be  completed  in  the  near  future,  this  problem  will 
become  even  more  critical  because  of  the  a oonmnanvi no'  lower  nressiires. 

—  —  y-  -  ^  O  A 


In  view  of  this  problem,  a  study  was  initiated  to  develop  a  flexible 
digital  computer  program  that  could  be  used  as  an  engineering  tool  to; 

1.  Study  the  effects  of  all  tubing  dimensions  upon  the 
system  response. 

2.  Determine  approximate  lag  times  for  use  in  determining 
more  realistic  test  schedules  and  minimum  stabilization 
times. 

3.  Determine  optimum  tubing  dimensions  within  given 
restraints. 

4.  Approximate  transducer  response  to  slowly  varying 
orifice  pressures,  including  complex  orifice  pressure¬ 
time  relationships. 

5.  Compare  relative  time  responses  of  different  tubing 
configurations . 

Several  theoretical  and  experimental  investigations  of  these  problems 
as  they  pertain  to  wind  tunnel  instrumentation  have  been  made.  Bauer 
(Ref.  1)  and  Kinslow  (Ref.  2)  have  made  detailed  studies  of  one-  and  two- 
tube  pressure  measuring  systems  of  constant  volume.  However,  since 
the  model  orifice  diameters  and  lengths  are  generally  fixed  by  the  test 
sponsor  and  the  transducers  are  installed  with  permanent  tubing  to  a 
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disconnect  panel,  a  three-tube  system  is  usually  employed  in  PWT.  The 
equations  of  Refs.  1  and  2  are  not  sensitive  to  the  order  of  the  tubes  in 
systems  in  which  more  than  two  tube  sizes  are  used.  The  theoretical 
results  of  Refs.  1  and  2,  therefore,  have  been  extended  to  include  three- 
tube  systems  in  such  a  manner  that  the  response  is  sensitive  not  only  to 
the  geometry,  but  also  the  order,  of  each  tube.  The  resulting  differ¬ 
ential  equation  relating  pressure  and  time  has  been  solved  for  both  step 
and  ramp  forcing  functions. 

The  purpose  of  this  report  is  to  present  the  theoretical  extension  and 
to  explain  the  utility  and  limitations  of  the  resulting  step  and  ramp  com¬ 
puter  programs.  Examples  are  given  to  serve  as  guides  in  the  employ¬ 
ment  of  these  engineering  tools.  Experimental  verification  of  second-tube 
diameter  optimization  is  also  included. 


SECTION  II 
SYSTEM  EQUATIONS 


The  computer  equations  are  based  upon  an  extension  of  the  theoretical 
analyses  of  unsteady,  compressible,  laminar  flow  in  cylindrical  tubing 
performed  by  Bauer  (Ref.  1)  and  Kinslow  (Ref.  2).  Bauer  considered  the 
distribution  of  volume  in  a  two-tube  constant  volume  system  but  neglected 
the  change  in  mass  flow  rate.  Kinslow  refined  the  theory  by  including  the 
effect  of  variation  of  mass  flow  in  a  two-tube  system.  To  account  for 
multidiameter  tubing  configurations,  Kinslow  suggests  replacement  of  all 
the  tubes,  with  the  exception  of  the  tube  next  to  the  transducer,  by  an 
equivalent  tube  having  the  same  volume,  mass  flow,  and  pressure  drop 
as  the  tubes  it  replaces.  The  equations  for  determining  the  geometry  of 
the  equivalent  tube,  however,  are  derived  assuming  constant  mass  flow 
in  the  replaced  tubes.  Hence,  the  equivalent  tube  is  not  sensitive  to  the 
order  of  the  replaced  tubes.  To  properly  study  the  flow,  and  in  particular 
to  optimize  the  diameter  of  any  given  tube  in  a  multidiameter  system,  the 
intuitive  conclusion  is  that  the  equations  describing  the  flow  must  be  sensi¬ 
tive  not  only  to  the  geometry,  but  also  to  the  order,  of  each  tube. 

The  analysis  is  based  upon  the  unsteady  form  of  the  Hagen -Poiseuille 
equation  of  laminar  flow  of  a  viscid  fluid  through  a  circular  tube  of  con¬ 
stant  diameter.  The  two-dimensional  flow  is  assumed  to  possess  axial 
symmetry,  with  the  velocity  distributed  parabolically  over  the  radius. 

The  following  assumptions  also  apply: 
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1.  The  process  is  isothermal. 

2.  Pressure  and  density  are  constant  over  any  cross 
section  of  a  tube. 

3.  Pressure  and  mass  flow  are  constant  across  any 
junction  of  two  tubes. 

4.  The  changes  in  pressure  and  density  along  a  tube 
are  linear. 

5.  The  times  and  lengths  required  for  complete  establish¬ 
ment  of  flow  are  small  compared  to  the  lag  times  and 
tubing  lengths  involved. 

6.  The  fluid  is  a  perfect  gas. 

7.  Acceleration  forces  are  negligible  compared  to  viscous 
forces. 

8.  Outgassing  and  leakage  rates  are  negligible. 

9.  Volume  of  tubing  transducer  system  is  constant. 

In  order  to  determine  the  time  response  of  the  transducer  pressure, 
a  differential  equation  relating  the  transducer  pressure  to  the  orifice  pres¬ 
sure  (forcing  function)  and  the  geometric  dimensions  of  the  system  is 
required.  The  methods  used  to  derive  this  equation  for  a  three-tube  sys¬ 
tem  are  similar  to  those  used  by  Kinslow  (Ref.  2)  in  deriving  the  two -tube 
equation.  The  pressure  measurement  system  for  which  the  equation  was 
developed  is  shown  schematically  in  Fig.  1  (Appendix  I).  A  detailed  deriva¬ 
tion  of  the  system  equations  is  included  in  Appendix  II.  The  system  equa¬ 
tions  are: 


DIFFERENTIAL  EQUATION 


pt2 


(1) 


GEOMETRIC  CHARACTERISTICS 


Km  =  5.16x10" 
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Kt  =  5.16  x  10' 


723  \  /  T, 


1.5 


^Tj  +198.6  /  y 524.4 
+  A2i2  (3B  -f  3DC)  +  A3  ^.a 


1 


IB  +  I/Xd, 


VHaa(3B  +  D 


(3DC)J  +  (^)(A)  |a2  4  (3C  +  1)  +  A3*3  (C2  +  3C) | 


K*  =  Km  .  10‘ 


(3) 

(4) 


STEP  SOLUTION 


lL 


K 


m 


2(Pj  +  5) 


loge 


2S(Pj  +  S)-gf 
e(2P  j  +  5) 


(5) 


RAMP  SOLUTION 


P 


Km  (4k) 


C,W  (a,  x)  —  W  (a,- 


(,C3W  (a,x)  +  W  (a,— x)  ; 


(6) 


SECTION  III 

APPROXIMATE  RAMP  SOLUTION 


The  ramp  solution  given  by  Eq.  (6)  is  quite  difficult  to  evaluate,  even 
when  using  a  digital  computer.  An  approximate  solution  that  is  quite  good 
for  a  great  many  cases  has  been  found.  Equation  (1)  can  be  rewritten  as 


The  third  term  on  the  right  of  Eq.  (7)  is  the  difference  between  the  two 
rates  of  pressure  change.  In  many  systems,  in  particular  a  Type  1  con¬ 
trol  system,  following  a  short  stabilization  period,  the  output  will  follow 
the  input  with  equal  rate  but  will  have  a  finite  difference  in  value  as  shown 
in  Sketch  2,  Appendix  II.  Hence,  the  third  term  on  the  right  can  be  neg¬ 
lected  in  cases  where  the  pressure  measurement  system-forcing  function 
combination  is  such  that  the  system  appears  as  a  Type  1  control  system. 
This  simplification  can  be  explained  in  another  way. 
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Equation  (7)  can  be  expanded  in  a  binomial  series  as 


kt\  ap.  kn 
p2  =  pa  _  v  (1  +  Ll - L  _  —5 

r  Nr  kJ  at  2 


1  +  -WM  — 


1  dP. 


Pt  dt 


Km3  (dP, 
2  \dt 


dP  dPt 
dt  dt 


Km/  1  I  / d2P  Km  d2P, 


Vdt 


Kt  dt' 


KJ  +  higher  order  terms 


(8) 


Hence,  neglecting  the  third  term  on  the  right  of  Eq.  (7)  is  equivalent  to 
neglecting  all  terms  of  order  2  and  higher  in  the  series  expansion  of  P^. 


The  approximate  ramp  solution  is  then  given  by 


r>2  _ 
Jr  — 


(9) 


Once  the  values  of  Km  and  Kq1  are  known  for  a  particular  system,  numer¬ 
ical  values  of  P  can  be  easily  hand  computed  using  Eq.  (9). 


SECTION  IV 

QUALIFYING  CONDITIONS 


The  applicability  of  Eq.  (1)  is  contingent  upon  the  validity  of  the 
assumptions  that  the  flow  is  laminar  and  the  acceleration  forces  are  small. 
From  Ref.  2,  the  Reynolds  number  condition  for  a  single  tube  is  given  by 


d3(Pt2~P2) 

64ft2£RTj 


<  2000 


(10) 


and  the  condition  for  neglecting  acceleration  is  given  by 


24  (•£/ d) 


Re  loge(P/Pt) 


»  1 


Expressing  the  junction  pressures  as 

p-  -  (itt)  (pi  +  S)  +  (btt)  p‘ 


=  (fn)<p.^(ft})p. 


(ID 


(12) 

(13) 
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the  six  qualifying  conditions  for  the  three -tube  system  are 

d,3 


Re,  = 


,pi +  8)1  - p’.} 


Re, 


d’3  /p. 

64/z2i2  RT,  \  12 

p’3/ 

d  3  { 

3  Jp2 

-  P’l 

64/ia-e3RT,  l  23 

/ 

24  (Vd.) 

al  - 

Re,  loge  /P,2  \ 

Vi + 

24  <t2/d2) 

a2  = 

Re2  lo8e 

24  (Vd,) 

a3  “ 

Re3  l°ge 

(14) 

(15) 

(16) 

(17) 

(18) 

(19) 


SECTION  V 

COMPUTER  PROGRAMS 


The  equation  lists  used  in  writing  the  step  and  ramp  computer  pro¬ 
grams  are  summarized,  for  convenience,  in  Appendixes  III  and  IV, 
respectively. 


5.1  STEP  PROGRAM 

The  step  program  is  primarily  used  to  determine  the  lag  time  follow¬ 
ing  a  step  change  in  orifice  pressure.  System  dimensions,  initial  pressure 
(initial  transducer  pressure  must  be  equal  to  initial  orifice  pressure), 
change  in  orifice  pressure,  and  desired  settling  error  are  inputs  to  the 
program.  The  standard  input  sheet  is  shown  in  Fig.  2.  The  time  lag  and 
qualifying  conditions  are  the  outputs.  The  step  program  has  been  written 
in  Fortran  II  for  use  in  the  PWT  Raytheon  520  computer. 

Since,  in  many  cases,  the  pressure  parameters  are  not  known,  a  com¬ 
parison  of  several  possible  tubing  systems  based  upon  geometry  and  fluid 
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properties  alone  is  desirable.  Applying  the  method  of  separation  of  vari¬ 
ables,  the  lag  time  can  be  written  as 

*L  =  Km  .  F(P)  (20) 

where  Km  is  a  function  of  geometry  and  fluid  properties,  whereas  F(p) 
is  a  function  of  the  pressure  conditions  alone.  Since  Km,  as  given  by 
Eq.  (2),  completely  characterizes  the  system,  comparison  of  different 
systems  can  be  accomplished  simply  by 'comparison  of  the  respective  Km 
values.  A  more  convenient  value  of  the  geometric  characteristic  is  K*, 
where 

K*  =  Km  .  10~3 
tL 

Generalized  plots  of  F(p)*  =  -  for  various  values  of  P,-,  6,  and  e  can  be 

K* 

generated  by  the  program.  The  actual  lag  time  can  then  be  determined 
from  these  plots  for  a  particular  system  (ith  system)  with  geometric  char¬ 
acteristic  K*  by 

tx  .  =  .  K*  (21) 

^  i  k*  1 

Typical  plots  of  t^/K*  are  shown  in  Fig.  3. 


If  the  diameter  and  length  of  the  first  and  third  tube  is  either  specified 
by  test  requirements  or  permanently  fixed,  as  is  usually  the  case  in  PWT, 
the  optimum  diameter  of  the  middle,  or  connecting,  tube  is  defined  as  the 
diameter  that  minimizes  the  value  of  K*.  The  smaller  the  tube,  the  larger 
the  flow  resistance  caused  by  viscosity  of  the  fluid.  As  the  tube  becomes 
larger,  the  volume  of  fluid  in  the  tube  that  must  be  moved  is  increased. 
Hence,  a  minimum  value  of  K*  always  exists.  The  computer  is  programmed 
to  optimize  the  second  tube  of  either  a  two-  or  three-tube  system  by  itera¬ 
tion  of  the  diameter  in  steps  of  0.  003  in.  The  output  contains  (d2)0pt> 

K*min,  and  the  K*  values  corresponding  to  d2  values  both  smaller  and 
larger  than  (d2)0pt  by  0.  003  and  0.  006  in. 


The  six  options  available  with  the  program  are  explained  in  Fig.  2. 
Typical  option  1,  5,  and  6  outputs  are  shown  in  Figs.  4,  5,  and  6, 
respe  ctively. 

The  viscosity  is  corrected  for  temperature  using  Sutherland's  formula 


P-i 


Mr 


/'I'r  +  1 98.6>\  /  tA% 
\T  +  198.6  /  \  T .  / 


(22) 
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where 


f- 

fJ-i  =  3.8  x  10  7  slugs/ft-sec 
Tr  =  524. 4°R 
Ti  =  (T  +  460)°R 
T  =  °F 


In  determining  the  qualifying  conditions,  the  initial  pressure,  Pj, 
is  used  in  Eqs.  (14)  through  (19),  resulting  in  worst  case  solutions. 


5.2  RAMP  PROGRAM 

Although  Eq.  (6)  is  an  exact  solution  and  could  be  programmed 
despite  the  relative  complexity  of  the  Weber  functions,  ft  was  apparent 
that  numerical  integration  of  the  differential  equation,  Eq.  (11-27),  is 
more  practical.  Since  the  initial  value  of  the  transducer  pressure  is 
known,  the  Runge-Kutta  method,  a  method  for  stepwise  solution  of  initial  - 
value  problems,  is  employed.  The  accuracy  of  the  Runge-Kutta  method 
is  dependent  upon  the  magnitude  of  the  increments  of  the  independent 
variable,  t.  The  program  is  set  up  so  that  the  initial  value  of  transducer 
pressure,  PQ,  may  differ  from  the  initial  value  of  orifice  pressure  P-p 
Hence,  if  Pf  =  Pi  (zero  slope)  and  PQ  $  Pi,  then  the  ramp  is  actually  a 
step  with  6  =  Pi  -  PQ.  To  determine  the  size  of  the  integration  incre¬ 
ments,  the  same  step  and  system  were  applied  to  both  programs  and  the 
increment  adjusted  until  the  solutions  coalesced,  as  shoyrn  in  Fig.  7. 

The  ramp  program  is  currently  written  in  Fortran  II  for  an  IBM 
7074.  The  standard  input  sheet  is  shown  in  Fig.  8.  Since  the  initial 
transducer  pressure,  PQ,  can  be  specified  independently  of  the  orifice 
pressure,  solutions  can  be  combined  by  using  the  transducer  pressure 
found  at  the  end  of  the  first  ramp  as  the  initial  transducer  pressure  for 
the  second  ramp  and  so  on.  In  this  manner,  complex  orifice  pressure¬ 
time  relationships  can  be  approximately  synthesized  using  a  series  of 
ramp  functions  and  the  transducer  response  determined. 

Although  5000  integrations  are  performed  per  period  t,  only  the 
results  of  integration  at  every  1/  10  r  are  printed  out.  A  typical  output 
appears  in  Fig.  9. 

In  determining  the  qualifying  conditions,  P  =  PQ  and  6  =  Pf  -  PQ 
are  used. 
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5.3  OTHER  PROGRAMS 

The  techniques  employed  for  the  ramp  forcing  function  may  be  applied 
in  like  manner  for  other  functions.  One  method  is,  of  course,  to  synthesize 
the  function  with  a  series  of  ramps  as  discussed  previously.  Or,  the 
Riccati  equation,  Eq.  (11-27),  with  the  corresponding  forcing  function  f(t) 
found  by  using  Eq.  (11-28),  may  be  integrated  numerically,  For  example, 
if  the  orifice  pressure  consisted  of  a  sinusoid  B  sin  lot  superimposed  upon 
a  steady  level  A,  the  forcing  function  would  be 

f(t)  =  —  j  (A  +  B  sin  cnt)2  -  KjcuB  cos  &>t[  (23 ) 


SECTION  VI 

DISCUSSION  OF  PROGRAM  UTILITY 


The  selection  of  the  better  one  of  the  two  programs  to  use  for  a  particu¬ 
lar  problem  depends  upon  what  forcing  function  is  to  be  used  and  what 
information  is  required.  Each  program  has  certain  unique  features.  The 
step  program,  for  example,  is  the  best  program  to  use  when  comparing 
systems  because  of  its  Option  5,  compute  K*  only.  In  addition,  it  must 
be  used  for  system  optimization.  The  ramp  program  achieves  its  greatest 
utility  when  used  to  synthesize  complex  forcing  functions.  Both  programs 
can  study  system  response  to  a  step  forcing  function  but  have  different  out¬ 
puts.  The  step  program  is  preferred  when  tj^  is  to  be  determined,  and  the 
ramp  program  is  preferred  when  a  plot  of  the  response  is  desired.  The 
remainder  of  this  section  is  devoted  to  the  discussion  of  typical  uses  of  the 
programs,  interpretation  of  data,  and  presentation  of  results. 


6.1  DIAMETER  OPTIMIZATION 

The  theoretical  lag  time  equations  have  been  shown  (Refs.  1  and  2)  to 
predict  measured  lag  times  to  within  10  percent  for  pressures  greater 
than  50  psf;  however,  in  the  optimization  of  a  two-tube  system  the  most 
important  consideration  is  not  the  accuracy  of  the  lag  time  but  rather  the 
accuracy  of  the  determination  of  the  diameter  of  the  second  tube  which  will 
minimize  the  lag  time.  Since  the  optimization  equations  were  not  verified 
in  either  Refs.  1  or  2,  a  laboratory  setup  was  employed  to  experimentally 
verify  optimization.  The  results  of  this  work,  given  in  Fig.  10,  show  the 
expected  error  in  K*  values  but  only  negligible  error  in  the  value  of  the 
optimum  diameter.  Figure  11  shows  some  of  the  measured  time  lag  data 
obtained  for  one  of  the  five  tubing  configurations  used  in  the  experimental 
verification. 


9 


A  t  DC-T  R  -66-225 


A  plot  of  the  optimum  diameters  and  K*  values  for  three  different 
orifice  diameters  is  shown  in  Fig.  12.  These  data  were  generated  by  the 
step  program  to  determine  orifice  size  required  to  limit  the  time  lag  to  a 
reasonable  value  and  to  determine  the  appropriate  d2  of  a  three -tube  sys¬ 
tem  to  be  used  during  an  actual  test  in  Tunnel  16T.  An  interesting  point 
is  demonstrated  in  Fig.  12.  As  shown,  for  a  tubing  system  in  which  the 
first  and  third  tubes  are  fixed,  the  diameter  of  the  second  tube  may  be 
smaller  or  larger  than  the  diameter  of  the  third  tube,  depending  upon  the 
other  tubing  dimensions.  This  is  contrary  to  the  generally  accepted  in¬ 
tuitive  deduction  that  purports  that  an  optimum  tubing  system  is  always  a 
stepwise,  uniformly  tapered  system,  i.  e.  ,  each  successive  tube  diameter 
is  larger  than  the  previous  diameter  as  one  goes  from  orifice  to  transducer. 

Since  the  use  of  an  optimum  diameter  tube  would  normally  involve 
complete  installation  of  new  tubing  for  each  test,  which  would  be  both 
costly  and  time  consuming,  the  question  arises  as  to  whether  one,  or  per¬ 
haps  two,  tube  diameters  might  represent  a  reasonable  compromise.  The 
step  program  can  be  used  to  logically  determine  a  best  compromise  as  the 
following  example  study  demonstrates. 

In  this  example  study,  model  orifice  diameters  of  0.  020,  0.  030,  0.  040, 
0.050,  0.  060,  and  0.  070  in.  were  assumed  representative  of  orifices  nor¬ 
mally  encountered  in  a  series  of  model  test  programs.  Each  diameter  was 
given  the  same  weighting,  i.  e.  ,  all  six  diameters  were  assumed  to  receive 
equal  usage.  The  remainder  of  the  two-tube  measurement  system  used  in 
the  study  is  explained  in  Fig.  13.  The  three  combinations  of  compromise 
diameters  studied  are  explained  in  Table  I. 


TABLE  I 

EXPLANATION  OF  COMBINATIONS  EMPLOYED  IN  COMPROMISE 
DIAMETER  STUDY 


Combination 

Compromise 

Orifice  Diameters  Compromised,  in. 

1 

A 

0.020,  0.030,  0.040,  0.050,  0.060, 

0.  070 

O 

O 

0.  020,  0.  030 

c 

0.  040,  0.  050,  0.  060,  0.  070 

3 

D 

0.  020,  0.  030,  0.  040 

E 

0.  050,  0.  060,  0.  070 
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The  deviation  from  the  optimum  case  for  a  particular  orifice 
diameter  (d^)  at  a  given  diameter  d.2  is 

AK*  =  K*  -  (K*)  opt  (24) 

and  the  total  deviation  is  given  by 

{<  -  «4>„P,}  (25) 

where  i  =  20  refers  to  d.^  =  0.  020  in.  The  compromise  diameter  for  a 
specific  compromise  combination  can  be  defined  as  the  diameter  that  re¬ 
sults  in  the  minimum  value  of  (AK*)^^  given  by 

|  <«*) 

The  most  desirable  combination  would  be  one  in  which  only  one  com¬ 
promise  diameter  is  used.  This  compromise,  for  the  example  study,  was 
compromise  A,  shown  in  Fig.  13,  with  a  compromise  tube  diameter-(d2) 
of  0,058  in.  and  a  total  deviation  of 

{<“*>  -  «7 

Referring  to  Fig.  13,  it  is  apparent  that  the  smaller  the  orifice 
diameter,  the  more  important  is  the  diameter  of  the  second  tube  because 
of  the  steepness  of  the  curve  and  larger  magnitudes  of  K*.  Combination  2, 
compromised  of  compromises  B  and  C,  resulted  in  two  compromise  diam¬ 
eters  -  one  to  be  used  for  the  smaller  orifices  (0.  020  and  0.  030  in.  ),  and 
the  other  for  the  larger  orifices.  These  diameters,  shown  in  Fig.  13,  are 
0.  048  and  0.  090  in.  ,  respectively.  The  total  deviation  resulting  from  the 
use  of  these  two  compromises  was 

{<“*>....!  }ml„  *  °'™ 

Combination  3  was  studied  to  determine  whether  the  0.  040 -in.  orifice 
should  be  grouped  with  the  0.020-  and  0.030  in.  orifices.  The  resulting 
diameters,  compromises  D  and  E  in  Fig.  13,  yielded  a  deviation  of 

{<“*> 

Thus,  combination  2  was  the  best  compromise  of  the  combination  studied. 


(K*)  -  CKj)  0f 


pt  l  mm 


(26) 
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6.2  RAMP  INPUT 


Transducer  responses  to  typical  ramps  of  positive  and  negative  slopes 
are  shown  in  Figs.  14  through  16.  In  each  case,  the  response  was  deter¬ 
mined  using  both  the  computer  and  the  approximate  equation.  As  shown 
in  Figs.  14  and  15,  the  approximate  equation  results  are  quite  close  to  the 
computer  results  whenever  the  input  and  response  slopes  are  nearly  equal. 
This  was  to  be  expected  since  the  approximate  solution  was  based  upon 
the  assumption  of  equal  slopes.  The  approximate  response  shown  in 
Fig.  16  is  indicative  of  the  discrepancies  that  can  occur  whenever  the 
approximation  assumptions  are  not  valid. 


Examples  of  synthesizing  complex  inputs  with  ramp  segments  are 
shown  in  Figs.  17  through  19.  The  first  (Fig.  17)  is  representative  of  a 
pressure  pumpdown  curve.  The  second  input  (Fig.  18)  is  typical  of  the 
pressure-time  relationship  resulting  from  traversing  a  nozzle  exit  at  a 
constant  rate.  Figure  19  is  an  example  of  synthesizing  a  periodic  wave¬ 
form.  In  this  example,  the  transducer  pressure  was  initially  set  at  the 
lowest  pressure  rather  than  the  average  value  of  the  input  pressure  in 
order  to  demonstrate  that  the  computer  program  would  give  both  the 


transient  and  steady-state  responses.  In  Fig. 
shift,  and  transient  are  quite  apparent. 


19,  the  attenuation,  phase 
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Fig.  1  Three-Tube  Pressure  Measurement  System 
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PWT 

7/12/65 


INPUT  INSTRUCTIONS 

THREE-TUBE  PRESSURE  LAG  COMPUTER  PROGRAM 
Step  Forcing  Function 


Item 


Instructions 


1 


2 

4 

5 


yXXXX  where:  y  =  1  for  letter  S 

y  =  2  "  "  T 

y  =  3  »  "  A 

y  =  4  "  "  B 

y  =  5  "  "  L 


Numeric 
Omit  "19" 
Option  Code  1 
2 

3 

4 

5 

6 


Error  =  psf  Ap 

ll  rt  nj  _  l  rr.  1  M 

-  7 o  r  ±Li<x  i  vaii 
"  =  psf 

"  =  %  Final  Val .  " 

K*  Computation  Only 
d2  Optimum  &  K*  Minimum 


psf 

psf 

%  Initial  Val 
%  Initial  Val 


8  FT3 

9-14  FT 

10&13  Zero  for  1-Tube 

11&14  Zero  for  2-Tube 

13  Zero  for  Option  6 

15  Unity  for  Options  5  &  6,  <  16 

16  Zero  "  "  "  " 

17  Unity  "  "  "  "  <  28 

18  Zero  "  "  "  " 


Fig.  2  Standard  Input  Sheet  for  Step  Program 
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INPUTS 

THREE -TUBE  PRESSURE  LAG  COMPUTER  PROGRAM 
STEP  FORCING  FUNCTION 


ITEM 


2.  Month 


3 .  Day 


4.  Year 


5.  Option  No. 


CONTACT _ 

SHEET _ OF 


7. 

Error 

8. 

Xducer  Volume 

9. 

Length  1 

10. 

Length  2 

11. 

Length  3 

12. 

Diameter  1 

13. 

Diameter  2 

15.  nA 


16.  n  values  of  Ap 


17 . nPi 


18. n  values  of  P 
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Fig.  3  Generalized  Plots  of  Program  Calculated  Pressure  Lag  Times -Step  Forcing  Function 
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40  60  80  100  200  400  600  1000  2000 


Initial  Pressure,  psf 


Fig,  3  Continued 
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6  8  10  20  30  40  60  80  100  200  300  500 

Initial  Pressure,  psf 


Fig.  3  Continued 
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.4 

.3 


40  50  60  80  100 


Initial  Pressure,  psf 


Fig.  3  Continued 


t  /K* ,  sec 
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Initial  Pressure,  psf 


Fig.  3  Continued 
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Initial  Pressure,  psf 


Fig.  3  Concluded 


10-18-65  PT0617 


3  -  Tube  Pressure  Lag 


CO 

tn 


Run  No. 
5.0000*00 

PI 

2.0000*03 

1.8000+03 

1.6000+03 

1.40001103 

1.2000+03 

1.0000+03 

8.0000+02 

6.0000+02 

4.0000+02 

2.0000+02 


Option 

1.0000+00 


tL 

1.007S+01 

1.1146+01 

1.2466+01 

1.4140+01 

1.6334+01 

1.9334+01 

2.3684+01 

3.0561+01 

4.3069+01 

7.2944+01 


K* 

9.1450+00 

Re  I 

9.6973+01 

8.7628+01 

7.8283+01 

6.8939+01 

5.9594+01 

5.0250+01 

4.0905+01 
3.1561+01 
2 .2216+01 
1.2872+01 


6 

1.0000+00 


Re  2 

4.9781+01 

4.4903+01 

4.0024+01 

3.5145+01 

3.0266+01 

2.5387+01 

2.0508+01 

1.5629+01 

1.0750+01 

5.8710+00 


s 

1.0000+02 


Re  3 

6.8063+01 

6.1308+01 

5.4553+01 

4.7798+01 

4.1043+01 

3.4288+01 

2.7533+01 

2.0778+01 

1.4023+01 

7.2675+00 


V 

3.5000-04 


al 

4.0141+04 
4.0141+04 
4.0140+04 
4.0139+04 
4.0138+04 
4.0136+04 
4.0132+04 
4.0125+04 
4. 010  8+04 
4.0037+04 


JL  I 

1.5000+01 


a2 

5.4013+05 

5.4010+05 

5.4011+05 

5.4013+05 

5.4012+05 

5.4012+05 

5.4009+05 

5.4008+05 

5.4002+05 

5.3979+05 


£2 

8.5000+01 


a3 

1.4699+05 

1.4699+05 

1.4699+05 

1.4698+05 

1.4698+05 

1.4698+05 

1.4697+05 

1.4696+05 

1.4692+05 

1.4674+05 


£  3 

3.4000+01 


<31 

3.9160-03 


<32  <33 

7.5000-03  5.4170-03 


Fig.  4  Typical  Step  Program  Output  -  Option  1 


AEDC-TR-66-225 


-  Tube  Pressure  Lag  IO-18-65  PT0617 


Run  No. 

Option 

K* 

V 

£  1 

£  2 

I  3 

dl 

d2 

d3 

1.0000+01 

5.0000+00 

9.1450+00 

3.5000-04 

1.5000+01 

8.5000+01 

3.4000+01 

3.9160-03 

7.5000-03 

5.4170-03 

2.0000+01 

5.0000+00 

9.3850+00 

3.5000-04 

1.5000+01 

8.5000+01 

4.0000+01 

3.9160-03 

7.5000-03 

5.4170-03 

3.0000+01 

5.0000+00 

9.71^9+00 

3 .5000-04 

1.6000+01 

8.5OOO+OI 

3.4000+01 

3.9160-03 

7.5000-03 

5.4170-03 

it.  0000+01 

5.0000+00 

2.9790+00 

1.3500-03 

1.5000+01 

1.000  oh  00 

0.0000+00 

3.9160-03 

3.9160-03 

7.5000-03 

5.0000+01 

5.0000+00 

1.1458+01 

1.3500-03 

1.5000+01 

8.5000+01 

1.0000+00 

3.9160-03 

7.5000-03 

3.9160-03 

Fig.  5  Typical  Step  Program  Output  -  Option  5 


3  -  Tube  Pressure  Lag  10-18-65  PT0617 


Run  No. 

Option 

(d2)opt 

K*min 

K®-6 

K*-3 

K*  +3 

K*  +  6 

V 

£  1 

£2 

£  3 

dl 

d3 

1.1000+01 

6.0000+00 

6.7765-03 

8.9689+00 

9.1276+00 

9.0127+00 

8.9843+00 

9.0499+00 

3.5000-04 

1.5000+01 

8.5000+01 

3.4000+01 

3.9160-03 

5.4170-03 

2.1000+01 

6.0000+00 

770949-03 

9.2951+00 

9.3729+00 

9.3040+00 

9.3369+00 

9.4220+00 

3.5000-04 

1.5000+01 

8.5000+01 

4.0000+01 

3.9160-03 

5.4170-03 

3.1000+01 

6.0000+00 

6.6880-03 

9-4573+00 

9.6177+00 

9.4991+00 

9.4795+00 

9.555&+00 

3.5000-04 

1.6000+01 

8.5000+01 

3.4000+01 

3.9160-03 

5.4170-03 

5.1000+01 

6.0000+00 

6.5955-03 

1.0883+01 

1.1009+01 

I.O895+OI 

1.0956+01 

1.1100+01 

1.3500-03 

1.5000+01 

8.5000+01 

1.0000+00 

3.9160-03 

3.9160-03 

Fig. 

6  Typica 

1  Step  Program  Output  —  Option  6 
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INPUTS 

THREE  -  TUBE  PRESSURE  LAG  COMPUTER  PROGRAM 
RAMP  FORCING  FUNCTION 

CONTACT _ 

ITEM  SHEET _ OF  __ 


1.  Project  No. 

2.  Month 

3 .  Day 

4.  Year 

5.  Run  No. 

6.  Xducer  Volume 

7.  Length  1 

8.  Length  2 

9.  Length  3 

10.  Diameter  1 

11.  Diameter  2 

12.  Diameter  3 

13.  Temperature 

14.  Initial  Orifice P. 

15.  Final  Orifice  P. 

16.  Initial  Xducer  P. 

17.  Time  Interval 

Item 


Instructions 


1 

2 

4 

6 

7-12 

8&11 

9&12 

13 

14-16 

17 


Example:  PT8002 

Numeric 

Omit  "19" 

it3 

ft 

Zero  for  1-tube 

Zero  for  2-tube 

°F 

psf 

sec 


Fig.  8  Standard  Input  Sheet  for  Ramp  Program 
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GO 

O 


PROJ.  NO. 

=  PT  -8002 

DATE  = 

6-16-65  RUN 

NO.  =  11 

A  1 

JL  2 

A  3 

d  1 

d2 

d  3  y 

2.0000+00 

3.5000+01 

3 .3000+01 

2.9200-03 

8.1700-03 

5.4200-03  3.5000-04 

Km 

KT 

K* 

T 

2.5856+03 

1.3336+03 

2.5856+00 

7.2000+01 

Re! 

Re2 

Re3 

a  i 

a  2 

a  3 

3.9154+01 

6.2276+00 

7.5186+00 

2.0274+03 

2.4322+05 

4.6452+04 

Po 

Pi 

Pf 

T 

1.0000+02 

1.0000+02 

2.0000+02 

1.0000+01 

TIPE 

P 

1. 0000+00 

9.5425+01 

2.0000+00 

9. 1997+01 

3 .0000+00 

8.9702+01 

4.00C0+00 

8.8533+01 

5 .0000+00 

8.8488+01 

6.0000+00 

8.9569+01 

7.0000+00 

9. 1773+01 

8 .0000+00 

9.5095+01 

9.Q0C0+00 

9.9524+01 

1 .0000+01 

1.0504+02 

Fig.  9  Typical  Ramp  Program  Output 
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dgj  in.  X  103 


Fig.  10  Comparison  of  Experimental  and  Theoretical  Values  of  K*  Showing  Optimization 
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4  I _ L 

40  50 


1000 


Initial  Pressure,  psf 


Fig.  1]  Comparison  of  Experimental  and  Theoretical  Time  Lags  Used  in 
Diameter  Optimization 
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0.0 


□  Approximate 


/ 

/ 


_i - 1 - 1 - 1 _ i _ 1  1  ■ _ i _ 

20  30  40  50  60  70  80  90  100 

Time  ,  sec 


Fig.  14  Response  to  Ramp  with  Positive  Slope 


Pressure 


Fig.  17  Response  to  a  Synthesized  Pumpdown 
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re  Probe  to  a  Synthesized  Nozzle  Pressure  Profile 


Pressure  ,  psf 
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APPENDIX  II 

DERIVATION  OF  SYSTEM  EQUATIONS 


BASIC  FLOW-PRESSURE  EQUATIONS 


The  basic  mass  flow  and  pressure  equations  required  for  the  develop¬ 
ment  of  the  system  equations  are  taken  directly  from  Eef.  2  without 
explanation.  If  the  flow  is  laminar  (Re  2000)  with  negligible  accelera¬ 
tion  forces,  mass  flow  is  related  to  pressure  by 


W 


dp 

dx 


P 


877-jaRT, 


(II- 1) 


For  a  single -tube  system,  the  orifice  pressure  and  transducer  pres¬ 
sure  are  related  by 


The  mass  flow  along  the  tube  is 


and  the  mass  flow  at  the  orifice  is  given  by 


W 


t 


dp  _A£  /dPt  + 
dt  2  \  dt 


(II- 2) 


(II- 3) 


(II-4) 


THREE-TUBE  EQUATION 


Since  the  third  tube  and  transducer  can  be  considered  as  a  single - 
tube  system  measuring  the  pressure  P23  at  the  orifice,  Eqs.  (II- 2)  and 
(II-4)  become 


(II-5) 

(II— 6) 


Applying  Eq.  (II -3)  to  the  second  tube  at  x  =  SL^*  wg3 


is  also  obtained  as 


w 


23 


dp 

dt 


12  -  K 


(II-7) 
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Setting  Eqs.  (II -6)  and  (II- 7)  equal,  solving  for  K  and  then  substituting 
this  expression  for  K  in  Eq.  (II-3)  with  P^.  =  Pig  and  P  =  P23>  yields  the 
mass  flow  in  tube  2  as 


w2 


V 


dp 

dt 


A2l2  /dP12 
~2  V  dt  + 


dP, 


dt 


A3i3  /dP23 
2  \  dt  + 


The  mass  flow  in  tube  2  can  also  be  found  from  Eq.  (II- 1)  as 


(II- 8) 


w 


2 


(II-  9) 


An  equation  for  the  square  of  the  junction  pressure  Pg3  can  now  be 
found  by: 


$ 

1.  Setting  Eqs.  ( II - 8 )  and  (II-9)  equal 

2.  Separating  variables 

3.  Integrating  from  p  =  Pig  at  x  =  0  to  p  =  Pg3  at  x  =  ig 

4.  Solving  for  P^gg 


The  result  is 

□  2 
JT  n  % 


Ja,V 

/dP  12 

dP23\ 

/  6 

\  dt 

dt  / 

dP, 


dt 


+  a2v 


dP23 

dt 


W, 


dP  A  J 


3*3*2 


/dP23  dp\j 
V  dt  dt  / 


+  P2 


2  dt  2 

/ 

Solving  for  Pg32  in  Eq.  (1-5)  and  setting  equal  to  Eq.  (II- 10)  yields 

p,»  - pl  -  K^r) tt +  (v+  Sf +  (nr +  nr)  "M 


(II- 10) 


A  i 

+  2K,  ((V  +  3  3 


A3 V)  d£  +  A,i,\  dp23 
3  /  dt  +  \  6  /  dt 


where 


K2  = 


8n[il2 


k3  = 


(II- 11) 


(II-  12) 


The  derivative  of  P23  must  now  be  eliminated.  Substituting  the  aver¬ 
age  values  of  8p/8x  for  tubes  2  and  3  into  Eq.  (II- 1),  the  tube  mass  flows 
are 


w2 


'  V  \ 

Qn(iRT1) 


(11-13) 
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-  m  ■  fe) 


(11-14) 


At  the  junction  of  tubes  2  and  3,  the  pressure  is  P23  and  W2  =  W3,  or 

/p«-p 


V  = 


Solving  for  P23 

where 

Differentiating  Eq.  (11-16) 

dP„ 


cp12  +  p 

23  C  +1 


c  - « 


1  /  dp.,  dp 

C  —  + 


dt  C  +  l  y  dt  dt  J 

dP23 

Substituting  for  — — —  in  Eq.  (II- 11)  and  collecting  terms  yields 

^  -  P2  -  Ka  ^  +  Kb  £ 


where 


Ka  = 


K, 


3(C  + 1) 


jAjij  (3C  +  1)  +  A343  (C2  +  3C)  j 


Kh  = 


2K, 


3lcT 


yp  |  A2i2  +  A3^(C2  +  3C  +  3)  +  3V(C  +  1)2| 


(II-15) 


(11-16) 


(II-17) 


(11-18) 


(11-19) 

(11-20) 

(11-21) 


Equation  (11-19)  is  the  differential  equation  relating  the  pressure  at 
the  transducer  to  the  pressure  at  the  junction  of  tubes  1  and  2.  The  inter¬ 
mediate  pressure  P23  and  its  derivative  have  been  eliminated. 

Elimination  of  pressure  P12  is  accomplished  in  basically  the  same 
manner  in  which  P23  was  eliminated.  Therefore,  the  steps  are  only 
briefly  outlined: 

1.  Find  w'x2  from  Eq.  (II- 8)  by  letting  x  =  0. 

2.  Find  W12  by  applying  Eq.  (II-3)  to  the  first  tube  at  x  =  SL\. 

Note  that  P^  =  Pt  and  P  =  P12  for  first  tube. 

3.  Equate  results  of  steps  1  and  2,  solve  for  K,  and  substitute 
the  resulting  expression  for  K  in  Eq.  (II-  3 )  applied  to  first 
tube  at  x  =  x,  resulting  in  an  equation  for  wi. 
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4. 

5. 

6. 

7. 


8. 


Find  W]_  by  applying  Eq.  (II- 1)  to  first  tube. 

Equate  Wi  expressions  of  steps  3  and  4,  separate  variables, 
and  integrate  from  p  =  patx  =  0top  =  pj2  at  x  =  SL\.  Solve 
for  P%i2- 

Solve  Eq.  (11-19)  for  P^^,  set  equal  to  P^i2  as  found  in 
step  5,  and  solve  for  P^2  -  p2.  The  pressure  Pj2  is  thus 
eliminated. 


Set  equation  for  wi  as  found  in  step  4,  with 


Sp 


equal  to  Eq.  (11-13)  at  p  =  Pi2*  Solve  for  P12  as  function  of 

dP 

Pt  and  P23.  Differentiate  P12  to  find  — .  Substitute 

dp.  dPia  dt 

expression  for  in  P^  -  P^  equation  of  step  6.  The 


dt 


derivative  of  Pi 2  is  thus  eliminated. 
dP 

The  derivative  — reappears  as  the  above  steps  are  per¬ 
formed  and  must  be  eliminated  from  the  P^  -  p2  equation 

dr 

resulting  in  step  7.  Express  — —  as  a  function  of  P-t  and  P  by 

dp  dt  dP 


using  Eq.  (11-18),  with 


replaced  by  the 


dt  ™-  dt 

found  in  step  7.  The  derivative  is  thus  eliminated 


expression 


The  resulting  equation,  containing  only  the  proper  independent  vari¬ 
able  (applied  orifice  pressure),  dependent  variable  (resulting  transducer 
pressure)  and  system  characteristics  (geometric  dimensions  and  fluid 
properties)  is 

P2  -  P2  =  KT  —  +  Km  4^-  (H-22) 

1  dt  111  dt 

where 


K.n 


87 Tfi 

T”  jvb 


1  \/t\r 

- — )  -tVHAA  (3B  +  1)  +  A2t2  (3B  +  3DC)  4  A3  i3  (3DC)] 

J  •  1  /  yA^  / 


+  1- 1— j  [Aa  t2  (3C  4  1)  4  A3t3  (C2  4  3C))J 


(11-23) 


K„ 


877 f-i  \(  1 


Ay)  [2A,  4  A2  i2  (3D  4  6)  4  A3t3  (3D  4  3B  4  6) 

4  1)  4  A323  (C2  4  3C)1 


3  (VB  +  i/VAj2, 

6V(B  +  1)]  ^)(-^)(tr)lA.M3C 


(£i)& 


2  ,  [A2  t2  4  A3  t3  (C2  4  3C  4  3)  4  3V(C  4  1)2]( 


(11-24) 
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(11-25) 


B  =  D(C  +  1) 


(11-26) 


Equation  (11-22)  is  identical  in  form  to  Eq.  (37)  of  Ref.  2.  The  geom¬ 
etric  characteristics  Kq1  and  Km  given  by  Eqs.  (11-23)  and  (11-24)  are 
different  in  that  they  now  contain  the  dimensions  of  a  third  tube.  For  a 
two-tube  system  (j(?3  =  0),  Eqs.  (11-23)  and  (11-24)  become  identical  to  the 
Kq  and  Km  equations  of  Ref.  2. 


SOLUTION  OF  DIFFERENTIAL  EQUATION 

Equation  (11-22)  can  be  rearranged  to  give 


dt  K„ 

wher'e  the  forcing  function  is  given  by 

P  2 

f(t)  =  ^  - 

K-  m 


-  =  f(t) 

a 

(11-27) 

KT  dPt 

(11-28) 

Km  dt 


Equation  (11-2  7)  is  a  nonlinear,  ordinary  differential  equation  of  first 
order  and  second  degree  with  constant  coefficients,  commonly  called  a 
Riccati  equation.  The  techniques  normally  employed  for  the  solution  of  this 
equation  depend  to  a  great  extent  upon  the  form  of  the  forcing  function  f(t). 


STEP  FUNCTION  SOLUTION 


For  a  step  in  pressure  forcing  function,  as  shown  in  Sketch  1, 


=  0  for  t  >  0  and 


f(t) 


Pt2  '  (Pj  +  5)2 
Km  Km 


(11-29) 


Pf  — _ 

_ d !_f _ 

Pressure  l 

r 

6 

„  L 

— 

\ 

\  ft 

V 

0  Time 


Sketch  1 
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The  Riccati  equation,  Eq.  (11-27),  becomes 

dp  P2  Pt2 

d  t  K  m  K-  m 


(11-30) 


If  a  particular  solution,  u,  is  known,  the  solution  of  this  form  of  the 
Riccati  equation  reduces  to  the  integration  of  a  linear  differential  equa¬ 
tion.  The  particular  solution  P  =  Pt  =  Pi  at  t  =  0  is  known. 


Letting 

P  =  P  +  I 

t  U 

(11-31) 

then 

1 

p-pt 

(11-32) 

and 

dP  _  1  du 

(11-33) 

dt  u2  dt 

t 

Substituting  Eqs.  (11-31)  and  (11-3  3)  in  Eq.  (11-30)  yields 


| 

1  du 

+  -L 

(p  -  'Y . 

(11-34) 

i 

u2  dt 

Km 

or 

du 

=  -L-j 

2Pt  u  +  1  } 

(11-35) 

dt 

Km< 

Separating  variables 


dt  =  Kn 


MM 

\2Pt  u  +  1/ 


(11-36) 


Now  at  t  =  0,  P  =  Pi  and  u  =  l/Pj-Pt,  and  at  t  =  t,  P  =  P  and  u  =  l/P-P^. 
Integrating 

1  # 


/*  =7_i?_ 

0  J  2pt  u  +  1 


(11-3  7) 


results  in  the  solution 


prpt 


(11-38) 


Equation  (11-38)  is  so  arranged  to  make  time  the  dependent  variable. 
Therefore,  the  lag  time,  or  the  time  required  for  the  transducer  pressure 
to  stabilize  to  within  a  specied  error,  e  ,  following  a  step  change  in  orifice 
pressure,  5,  can  be  expressed  as 


Km 

2(P;  +  <5) 


l°ge 


25(Pi  +  §)  -  8  ( 
e(2Pj  +  5) 


(11-39) 
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RAMP  FUNCTION  SOLUTION 


For  a  ramp  in  pressure  forcing  function,  as  shown  in  Sketch  2, 

— L  -  _J - L  for  t  >  0  and  f(t)  =  kt2  +  mt  +  n  (II~40) 

dt  t 


where 


(11-41) 


(11-42) 

(11-43) 


Pressure 


Through  the  use  of  a  transformation,  the  first  order  nonlinear  Riccati 
equation  can  be  transformed  into  a  second  order  linear  differential  equa¬ 
tion  that  can  be  solved.  Such  a  transformation  is 

Y  =  e  / Ps(t)dt  (11-44) 

where  g(t)  is  a  function  of  the  independent  variable,  time,  that  will  make 
the  differential  equation  linear. 

Differentiating 

i>L  =  e  /pe  dt  —  [  f  pg  dt]  (H-45) 

dt  dt  ’ 

or 

-J-  =  Pgy  (11-46) 

or 

P  =  -i-  if  (11-47) 

yg  dt 
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Differentiating  Eq.  (11-47) 


dp  =  1 

dt  g 


J_  _d_  /l_  dy\  _  J_  /J_  dy\  dg 
g  dt  \y  dt/  g2  \y  dt/  dt 


(11-48) 


Since 


d  /  1  dy\  _  1  d2y  1  dy  dy 

dt  \y  dt  /  y  dt2  y2  dt  dt 

=  L  ."Lx  _  /L  dy\2 

y  dt2  \y  dt  / 


Eq.  (11-48)  becomes 


dp  _  1  d2y  _  _  /  dy\  _ 1_  dy  dg 

dt  gy  dt2  gy2  \dt  /  g2y  dt  dt 


Substituting  Eqs.  (11-47)  and  (11-50)  into  Eq.  (11-27)  results  in 


1  d2y 


_ L_  /ilV _ L_ 

gy2  \dt  /  g2y 


dy  dg  _ 

dt  dt  K 


—  (—)  = 
ng2y2  \dt  / 


=  f(t) 


(11-49) 


(11-50) 


(11-51) 


If  g(t)  is  chosen  to  be  equal  to  1/Kmj  which  makes  =  0,  Eq.  (11-51) 

dt 

reduces  to  the  linear  equation 


d2y  f  (t)  y 


and  the  inverse  transform  is 


p  =  dy^ 

y  dt 


(11-52) 


(11-53) 


If  the  mt  term  of  the  f(t)  function  is  eliminated,  Eq.  (5  2)  will  be  in  the 
form  of  a  differential  equation  for  which  a  solution  is  known.  Therefore, 
let 


(4k)%  2k 


therefore 


Using  the  transformation,  f(t)  becomes 

f(t)  =  f(x)  =  ^5-  -  —  +  n 
2  4k 

Substitution  of  Eqs.  (11-55)  and  (11-56)  into  Eq.  (11-52)  results  in 


-  a]  y  =  0 
L4Km  J 


where 


m*-  4kn 
8Kmk3/2 


(11-54) 


(11-55) 


(11-56) 


(11-57) 


(11-58) 
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Equation  (11-57)  is  a  form  of  Weber's  equation.  Its  solutions  are 
parabolic  cylinder  or  Weber-Hermite  functions  which  can  be  found  tabu¬ 
lated  in  Ref.  3  or  found  by  the  methods  of  Ref.  4.  The  solution  is 

y  =  C,  W(a,x)  +  C2  W(a,  -x)  (11-59) 


where  the  Weber  function  is  given  by 


where 


(11-60) 

(11-61) 


3/2)  fr 


+  (a3  -  13/2a)  + 


(11-62) 


G,  =  ir  (1/4  +  i  a/2)  I 
G.  =  T  (3/4  +  i  a/2)  | 

Therefore 

P  =  Kra  (4k)  ^ 

where 

W'(a,x)  =  llWi(a,^) 
dx 

No  loss  in  generality  occurs  if  C2  =  1,  therefore 

P  =  Km  (4k) ^ 


C jW  ,  (a,  x)  -  W'(a,-x)-\ 
CtW  (a,  x)  +  W(a,-x)  j 


C,W  '  (a,  x)  -C2W"(a,-x)1 
C,W  (a,  x)  +  C2W  (a,  — x)  j 


(11-63) 

(11-64) 

(11-65) 

(11-66) 

(11-67) 
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APPENDIX  III 

STEP  FORCING  FUNCTION  EQUATION  LIST 


AREAS 

A, 

A2 

a3 


DIMENSIONLESS  CONSTANTS 


C 


n 


/v\  (<\2 
WW 


0.7854  dj2 

0.7854  d22 
0.7854  d32 


C  =  0  if  i3  =  0 


D  =  0  if  t2  =  0 


B  =  I)(C  +  1) 


JUNCTION  PRESSURES 

P‘-  =(5TT)(Pi  +  5)+(lTT)P. 

P..  -  (to)  (Pi -  s)  *  (w)  p. 


REYNOLDS  NUMBER 


(III- 1) 
(III- 2) 
(III-3) 


(III-4) 

(HI-  5 ) 

(III- 6) 


(III- 7) 
(III- 8) 


(III- 9) 
(III- 10) 

(III-ll) 
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ACCELERATION  CONDITIONS 


•i 


24  (Vdi> 


Bej  logc  p, 


P.  +  S 

24(l2/d2) 

|  Re2  l°Se  P23 

P12 

24(Vd,) 
Re3  logc  Pj 


GEOMETRIC  CHARACTERISTICS 


Km  =  5.1§  x  TO-6  1-5  (-*— )  Pv)  Ua,  ^ 

f  \Tj+ 198.6/  \524.4/  \B  +  l/  \  d.V 


+  A2 /2  (3D  +  6)  +  A3  i3  (3D  +  3B  +  6)  +  6V(B  +  l)f 


AJ2  +  A.*,  (3C  +  1)  +  A ,L  (C2  +  3C) 


A2 t2  +  A3  lz  (C2  +  3C  +  3)  +  3V(C  +  l)2 


K*  =  Km  .  10' 


LAG  TIME 


lL  = 


K, 


2(P  .  +  8) 


!°ge 


2S(Pi  +  S)~5f 


e  (2F+S) 


UNIT  CONVERSIONS 


Test:  If  P.  +  8  <  0 


then  tL  =  0 


£  (Percent  Final  Value)  _  p  Options  2  and  4 
100  1 


^  S  (Percent  Initial  Value)  ^  „  .  , 

o  =  - -  .  Fj  Options  3  and  4 


(III- 12) 


(III- 13) 


(III- 14) 


(III- 15) 


(III- 16) 


(III- 17) 
(III- 18) 
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OPTIMIZATION  STARTING  VALUE  (OPTION  6) 


TEMPERATURE 


T,  =  T  +  460 


(III- 19) 

(III- 20) 


(III-21) 
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APPENDIX  IV 

RAMP  FORCING  FUNCTION  EQUATION  LIST 


AREAS 

A,  =  0.7854  d,2 
A2  =  0.7854  d22 
A3  =  0.7854  d32 

DIMENSIONLESS  CONSTANTS 

c'(t)  (t)’ c  - 0 - 0 

«  (h\  /Ai\2  „  .....  _ 

u  JJ  =  0,t**  =  u 

B  -  D(C  +  1) 

JUNCTION  PRESSURES 

p“  -  (sfr)  <pi  *  -  po> + 
p-  -  (stt)  <pi +  p<  -  p.> +  (frf)  p< 

REYNOLDS  NUMBER 


Re,  = 


1.1905  x  io5  (r~){(pi +  pf  -  po)2  - p-2} 


Re2  = 


Re.  = 


1.1905  x  10s 


1.1905  x  10s 


© 

)  (P..2 

\i3  /  '  “ 


(P122  -  P232) 


P  2) 


(IV- 1) 
(IV-2) 
(IV-3) 

(IV-4) 

(IV- 5) 
(IV-6) 


(IV- 7) 
(IV- 8) 


(IV-9) 
(IV- 10) 
(IV- 11) 
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ACCELERATION  CONDITIONS 


TEMPERATURE 


24  (Vd.) 

P12  \ 

;+pf-po/ 

24  (Vd2) 

Re2  loge  /PjA 

\pJ 

24  (l3/d3) 

Re3  /  P0\ 

\P  23  / 


Rei  log, 


6 


T,  =  T  +  460 


(IV- 12) 


(IV  13) 


(IV- 14) 


(IV- 15) 


GEOMETRIC  CHARACTERISTICS 


Km  =  5.16  x  10 )  AA  '  p-)  (-^)  j2At  K 
V^  +  198.6/ V524.4/  |_\B  +  l/\d14/  ( 

+  A2U2  (3D  +6)  +  A3^3  (3D  +  3B  +  6)  +  6V(B  +  1)| 


V,  +  A3I3  (C2+3C  +  3)  +  3V(C  +  1)2;: 


(IV- 16) 


KT  =  5.16  xlO-6  (  ■  7-2--- )  (-^A  '  (^A  (A)  |a3  (3B  +  1) 

VTj  +  198.6  /  '524. 4'  [\B  +  1/  'df/  < 

+  A2  2  (3B  +  3DQ  +  A3  i,  (3DC)  j  +  (-^j  (-^j  j  A2  Jt  (3C  +  1) 


+  A3  ia  (C2  +  3C)J 


(IV- 17) 


K*  =  Km  .  10“3 


(IV -18) 
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FORCING  FUNCTION  COEFFICIENTS 


■A 


DIFFERENTIAL  EQUATION 


dP 

dt 


— —  -=  kt2  +  mt  +  n 

^  m 


INTEGRATION 


Initial  Conditions:  P  =  P  at  t  =  0 

o 

dP/dt  =  0  at  t  =  0  if  PQ  =  P. 
Integration  Increment:  At  =  0.0002  r 
Printout  Increment:  t  =  O.lr,  0.2r,  .  .  .  l.Or 


(IV- 19) 
(IV-20) 
(IV-21) 


(IV-22) 

(IV- 23) 

(IV-24) 

(IV-25) 

(IV-26) 


* 


! 


.  ; 


. 


! 
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